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Abstract 

In this paper we propose and analyze a Discontinuous Galerkin 
method for a linear parabolic problem with dynamic boundary con¬ 
ditions. We present the formulation and prove stability and optimal 
a priori error estimates for the fully discrete scheme. More precisely, 
using polynomials of degree p > 1 on meshes with granularity h along 
with a backward Euler time-stepping scheme with time-step At, we 
prove that the fully-discrete solution is bounded by the data and it 
converges, in a suitable (mesh-dependent) energy norm, to the exact 
solution with optimal order h p + At. The sharpness of the theoretical 
estimates are verified through several numerical experiments. 


1 Introduction 

In this paper we present and analyze a Discontinuous Galerkin (DG) method 
for the following linear parabolic problem supplemented with dynamic bound- 
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ary conditions on Tj: 

dtu = Au + /, in D, 0 < t < T, 

d n u = — cm + (3Ayu — Xdtu + g, on Id, 0 < t < T, . . 

periodic boundary conditions, on Id, 0 < t < T, 

„ u \t=o = Uo, in fi. 

Here the domain H and the subsets Tj C dQ, i = 1,2, are depicted in 
Figure[lJ Ap is the Laplace-Beltrami operator, d n u denotes the outer normal 
derivative of u on Id, g is a given function and a,/3,X are suitable non¬ 
negative constants. 

Dynamic boundary conditions have been recently considered by physi¬ 
cists to model the fluid interactions with the domain’s walls (see, e.g., 
he \m mg). Despite the practical relevance of this kind of boundary con¬ 
ditions from a modeling point of view and the intense research activity to 
understand their analytical properties, see, e.g., m mm, the study of 
suitable numerical methods for their discretization is still in its infancy. To 
the best of our knowledge, the only work along this direction is [5], where 
the authors analyze a conforming finite element method for the approxima¬ 
tion of the Cahn-Hilliard equation supplemented with dynamic boundary 
conditions. Motivated by the flexibility and versatility of DG methods, here 
we propose and analyze a DG method combined with a backward Euler time 
advancing scheme for the discretization of a linear parabolic problem with 
dynamic boundary conditions. The main goal of the present work is the 
numerical treatment of dynamic boundary conditions within the DG frame¬ 
work. Here we consider just a linear equation. However, our results aim to 
be a key step towards the extension to (non-linear) partial differential equa¬ 
tions with dynamic boundary conditions, as, for example, the Cahn-Hilliard 
equation. In this context, we mention DG methods have been already proved 
to be an effective discretization strategy for the Cahn-Hilliard equation as 
shown in PS where the authors constructed and analyzed a DG method 
coupled with a backward Euler time-stepping scheme for a Cahn-Hilliard 
equation in two-dimensions, cf. also M- 

The origins of DG methods can be backtracked to El EE where they 
have introduced for the discretization of the neutron transport equation. 
Since that time, DG methods for the numerical solution of partial differen¬ 
tial equations have enjoyed a great development, see the monographs mm 
for an overview, and |3j for a unified analysis of DG methods for elliptic 
problems. In the context of parabolic equations, DG methods in primal 
form combined with backward Euler and Crank-Nicholson time advancing 
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techniques have been firstly analyzed in mm, respectively. DG in time 
methods have also been studied for parabolic partial differential equations, 
see, for example, US 0 ES [ 20 ] and the reference therein; cf. also [28, [27] for 
the hp-ve rsion of the DG time-stepping method. 

The paper is organized as follows. In Section [2] we introduce some useful 
notation and the functional setting. Section [3] is devoted to the introduction 
and analysis of a DG method for a suitable auxiliary (stationary) problem. 
These results will be then employed in Section [4] to design a DG scheme to 
approximate the linear parabolic problem with boundary conditions and to 
obtain optimal a priori error estimates for the fully discrete scheme. Finally, 
in Section [6] we numerically assess the validity of our theoretical analysis. 

2 Notation and functional setting 

In this section we introduce some notation and the functional setting. 

Let D C M 2 be an open, bounded, polygonal domain with boundary 
T = dD. On D we define the standard Sobolev space H S (D), s = 0,1, 2,... 
(for s = 0 we write L 2 (D) instead of H°(D)) and endow it with the usual 
inner scalar product (■,-).ff s (D)) and its induced norm ||-|| #«(£>), cf. [I]. We 

also need the seminorm defined by • I „.<D) = (£ w =, I|9 "(-)IIl»(0)) 1/2 . 

We next introduce, on T, the Laplace-Beltrami operator. We first define 
the projection matrix P = I — n ® n = (dij — rii'nj)f ]=l . where n is the 
outward unit normal to D, a (gib = (aibj)ij is the dyadic product, and Sj-j is 
the Kroneker delta. We define the tangential gradient of a (regular enough) 
scalar function u : T —> M as Vpu = PVit. The tangential divergence of 
a vector-valued function A : T —> M 2 is defined as divr(A) = Tr((VA)P), 
being Tr(-) the trace operator. With the above notation, we define the 
Laplace-Beltrami operator as Apu = divp(Vru). 

We next introduce the following Sobolev surface space 

H S (T) = {ve H S ~\T ) | V r u € [iL s_1 (T)] 2 }, s > 1, 

cf. [2, with the convention that H°(T ) = L 2 (T), L 2 (Y) being the standard 
Sobolev space of square integrable functions (equipped with the usual inner 
scalar product (-,-)r and the usual induced norm 11• 11z> 2 (r))- W e equipped 
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the space H S (T) with the following surface seminorm and norm 

M/r»(r) = ||V r u||fr—i(r) Vn € H S (T), s > 1, 

IMI-ff s (r) i(r) l^l_ff s ( r ) ^ ^ (r)j s > 1, 

respectively. In m Lemma 2.4] is proved that the above norm is equivalent 
to the usual surface norm present in literature [122] . which is defined in local 
coordinates after a truncation by a partition of unity. 

Next, for a positive constant A, we introduce the space 

H a x{D,T) = {v€ H S (D) : Au, r G H s ( L)}, a > 0, 
and endow it with the norm 


W))' 

As before, for s = 0 we will write HUD,T) instead of H®(D,F). Moreover, 
to ease the notation, when A = 1, we will omit the subscript. 

Finally, throughout the paper, we will write x < y to signify x < Cy , 
where C is a generic positive constant whose value, possibly different at any 
occurrence, does not depend on the discretization parameters. 

3 The stationary problem and its DG discretiza¬ 
tion 

Let O = (a, b) x (c,d) C R 2 be a rectangular domain and let ri,r 2 be the 
union of the top and bottom/left and right edges, respectively, cf. Figure [I] 
We consider the following Laplace problem with generalized Robin boundary 
conditions: 

—A u = /, in O, 

d n u = —olu + (3Aru + g, on Fi, (2) 

periodic boundary conditions, on F' 2 , 

where a,/3 are positive constants, and / G L 2 (0), g G L 2 (Ti) are given 
functions. 

Defining the bilinear form a(u,v) : IL 1 (0, Ti) x iL 1 (0, Ti) —> M as 
a(u, v) = (Vu, Vu) L 2 (Q) + /3(V r u, Vru) L 2 (ri) + a(u, v) L 2 {Tl) , 


u 


h x(d,t) 


|2 
I H ! 


(D) 


+ ^ll u |I 
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( 3 ) 


the weak formulation of Q reads: find u 6 i^ 1 (f7,ri) such that 

a{u,v) = {f,v) L 2( n) + (g,v) L 2 (ri) Vu E i). 

The following result shows that formulation ([3]) is well posed. 

Theorem 3.1. Problem |$|] admits a unique solution u E H 2 (Q, Id) satis¬ 
fying the following stability bound: 

IMI// 2 (n,ri) < II/IIl 2 (Q) + Hs'IIlRi’i)- (4) 

Moreover, if f E H s ~ 2 (Pl) and g G H s ~ 2 (Ti), s > 2, then u E H S (Q,T i) 
and 

IMI.ff s (f2,ri) ^ ll/llj^-Rn) + Ibllip-Rri)- (5) 

Proof. The existence and uniqueness of the solution are proved in [17, The¬ 
orem 3.2]. The proof of the regularity results is shown in [IT] Theorem 
3.3-3.4], The same arguments used in [l?f Theorem 3.3-3.4] apply also in 
our case thanks to periodic conditions. □ 

Remark 3.2. We observe that the forthcoming analysis holds in more general¬ 
shaped domains and/or more general type of boundary conditions provided 
that the exact solution of the differential problem analogous to 0 satisfies 
a stability bound of the form of Q). 

3.1 Discontinuous Galerkin space discretization 

In this Section we present a discontinuous Galerkin (DG) approximation of 
problem <§• 

Let Th be a quasi-uniform partition of 17 into disjoint open triangles T such 
that Q = UWe set h = max{diam(T), T G Th}- For s > 0, we define 
the following broken space 

H s (Th ) =K L 2 (n) : V\ T G H S (T, dT), T G T h \, 

where, as before, H°(Th) = L 2 (Th). For an integer p > 1, we also define the 
finite dimensional space 

V p (T h ) ={«G L 2 (n) : v ]T G P p (T), T G Th} C H s (%), 

for any s > 0. An interior edge e is defined as the non-empty intersection 
of the closure of two neighboring elements, i.e., e = T\ n T%, for T\, G Th- 
We collect all the interior edges in the set fP. Recalling that on f 2 C dQ we 
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impose periodic boundary conditions, we decompose T 2 as T 2 = Tj" UT^, 
cf. Figure [l] (left), and identify with T^", cf. Figure [I] (right). Then 
we define the set £^ 2 of the periodic boundary edges as follows. An edge 
e G S^ 2 if e = dT n 8T + , where G Th such that dT± C rf , cf. Figure [l] 
(right). We also define a boundary edge er 1 as the non-empty intersection 
between the closure of an element in Th and T i and the set of those edges 
by Sfr 1 - Finally, we define a boundary ridge r as the subset of the mesh 
vertexes that lie on Ti, and collect all the ridges r in the set Clearly, 

the corner ridges have to be identified according to the periodic boundary 
conditions (cf. Figure JTJ right). The set of all edge will be denoted by £h , 
i. e .,4=£!u£['u£f. 



Figure 1: Example of a domain and an admissible triangulation Th. (left). 
On the right, we highlight the edges e G £^ 2 with red lines. 


For v G H s (Th), s > 1, we define 

r6Th 

Next, for each e G £^ 2 we define the jumps and the averages of v G 
II'(Ti.) as 

He = (w + )n^ + (n _ )n“ and {v} e = ^(n + + n“), 

where = v^± and is the unit normal vector to e pointing outward of 
T±. For each e G £^ r we define 

[v\ e = v\ e n e , {v} e = v\ e , v g H l (Th)- 

Analogously, for each r G 7 ^ , we set 

Hr = (u + (?’))n+ + (n“(r))n“ and {w} r = ^(u + (r) + v~(r)), 
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where, denoting by eF the two edges sharing the ridge r, v±(r) = v\ e ±(r) 
and is the unit tangent vector to Id on r pointing outward of eA. The 
above definitions can be immediately extended to a (regular enough) vector¬ 
valued function, cf. [3j. To simplify the notation, when the meaning will be 
clear from the context, we remove the subscripts from the jump and average 
operators. Adopting the convention that 

(v,w) £h = ^2(v,w) L2{e) , (^,v) 7l r 1 = ^2 £(r)r]{r) 

e££h re^ 1 

for regular enough functions v, w, £, 77 , we introduce the following bilinear 
forms 

Bh(v,w) = Y (Vu, Vw) T - (M,{Vw}) £ o - (H,{Vu}) £ o +<t([u], [w])go 

T&T h 

- (H, {Vtn}) r 2 - (H, {Vu}) r 2 + a([v], [w]) r 2 

°h °h °h 

and 


b h (v, w) = (V r u, V r w) n -([u], {V r w}) n -(H, {V r u}) 17 +ct([u], [w]) r x , 

°h K'h h ’^h 

for all v,w G H 2 (Th )• Here a = ?, being 7 a positive constant at our 
disposal. We then set 


A h (u, v) = B h (u, v)+a (u, v) L 2 ( Tl) + (3 b h (u, v). ( 6 ) 

The discontinuous Galerkin approximation of problem ([ 2 ]) reads: find 
Uh G V p (Th ) such that 


A h (u h , v h ) = (/, v h ) L 2 ( U) + ( g , Ufc) L 2 ( Fl) \/v h G V p (%). (7) 

In the following we show that the bilinear form Ah(-, •) is continuous and 
coercive in a suitable (mesh-dependent) energy norm. To this aim, for w G 
H s (Th), we define the seminorm 


\w 


III 2 

HI B h 


W \HHT h ) 


+ MMII 


L2( £ o u£ r 2) 


+ -|I{Vuj}|| 
a 


£2 K^l 2 ) 


and the norm 


\w L = 


M\l h +a\\w " 2 


L 2 F 1 ) 


+ /3|w| ^(< 1 ) +/3c7|l[u;]ll i 2 (< 1 ) 


P 


a 


ll{VrMlli 2(<1) , ( § ) 
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where we adopted the notation 

IHli 2 (^) = 5ZlHli 2 (e)’ IMl' 2{ ^r 1} = J2 IMIW)- 

eE£: ' 1 eS^ 1 

Reasoning as in [2J, it is easy to prove the following result. 

Lemma 3.3. It holds 


■Ah{v,w) < IIMIUIMII* Vv,w £ H 2 (Th). (9) 

Moreover, for 7 large enough, it holds 

llklll" <Mv,v) Vu G V p (%)- ( 10 ) 

Proof. Let us first prove ©• The term Bh(-,-) can be bounded by Cauchy- 
Schwarz inequality as in [ 2 ]. Also the term bh(-, •) can be handled using the 
Cauchy-Schwarz inequality: 


\bh(v,w )| = (V r u,V r w),r 1 - (k],{V r w}) 17 

°h '^h 

-(H,{V r u}) r, +o-([u], [in]) r, 

+ -H{ VrU }lli2 (7e ri)) X 


< 

r^j 


1 \ 1/2 

w|^ (£ r 1 , + .||W|| 2 I3(K r 1) + -||{Vrw}||^ (R r 1) N 


and ([9]) follows employing the definition ([ 8 ]) of the norm ||| • H^. 

We now prove (10). As before the term Bh(-,-) can be bounded as in [ 2 ]: 
using the classical polynomial inverse inequality [ 6 ] we obtain 


I III 2 I |2 
\V I- < l/ ” 1 


cr|| |^| ||7 - T r 2 <B h (v,v) 


+U|| L“J"L2 {£ 0 U£ X2 ) 

for all v £ V p (Th )• The term 6 /,(-, •) can be estimated as follows: 


bh(v,v) > - 2 


|2 


(M,{Vrv}) w ri 


+ a||H 


'lhK 1 )' 


Employing the arithmetic-geometric inequality we get: 

(M,{V r u}) 7 ? r 1 < lk 1 / 2 k]|| r 2 f 7 ? r 1 J|{cj~ 1 / 2 Vru}|| 


< —cr|| [u] 


v \\\ L \n^)M a ' Vr ' u lllLR^) 
l' 2(< i )+ 4 6a- 1 ||{ v ru}|| 2 2 (.), 


8 









for a positive e > 0 . 
inverse inequality 


Finally, estimate (10) follows using the polynomial 


h ||{VrMH 


2 

lhK 1 ) 


< 


v 


2 


and choosing 7 sufficiently large. 


W € V p (T h ) 


□ 


The following result shows that problem 0 admits a unique solution and 
that the Galerkin orthogonality property is satisfied. The proof is straight¬ 
forward and we omit it for sake of brevity. 


Lemma 3.4. Assume that 7 is sufficiently large. Then, the discrete solution 
Uh of problem 0 exists and is unique. Moreover, formulation 0 is strongly 
consistent, i.e., 

A h (u-u h ,v) = 0 VveV p (T h ). (11) 

For v E H S (Q,T 1 ), s > 2, let ijfv be the piecewise Lagrangian interpolant 

of order p of u on Th- Note that (IpU) |p 1 interpolates u on the set of degrees 

r 

of freedom that lie on Sffi. By standard approximation results we get the 
following interpolation estimate. 


Lemma 3.5. For all v E H s (Q,Ti), s> 2, it holds 

Proof. Using the definition ^ of ||| • |||^ norm and that IpV(r ) = v(r ) for all 
r E , we get 


\\v - IpV\\f = HI?; - IpV\\Y 
11 p in* in p \»g h 


+ a\\v - ffiv || L 2 ( ri ) + f3\v - Ipvl^^y ( 12 ) 


Expanding the first term at right-hand side and using the multiplicative 
trace inequalities 


,2 <fc-Hl-.ll2 , U..I2 


1 L 2 (e h ) 

2 


I L 2 (Q) 


+ h\v\ 


2 


II^IIl 2 ^) ~ h Miji(n) + h\ v \ h 2 (Q)’ 


cf. [25] , we get 


lll v - = l u - T P V \HHC1) + ♦ - 

+ ^|{V(,-I»}||^ 2( ^ 2) 

<h 2 ||?; - IpV || 2 2 ( n) + \ v ~ ^p v lffi(n) + h2 \ v ~ -^p' u Ih 2 (q)- 
Using standard interpolation estimates [23] we get the thesis. □ 
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Now we show that the discrete solution Uh of 0 converges to the weak 
solution of <©■ 

Theorem 3.6. Let u £ H (^”^5 ri)^ s ^ 2, be the solutzon of the problem 
and let Uh be the solution of the problem a- Then, 

||«-«/»||i a (n,r 1 ) + h|||“-“h|||* < / l mm( *jH-i)|| lt || H , (n) r i ) j 

provided 7 is chosen sufficiently large. 

Proof. By the triangular inequality we have 

111“ - “hill* < 111“ - ^“lll* + 111-#“ - lihlll^. 

We first bound the second term on the right-hand side. Combining the 
Galerkin orthogonality 0 with the continuity and the coervicity estimates 
(j9j-(10), we obtain: 

III IpU - u h \\\l < A h (IpU - u h , IpU - u h ) 

= A h (IpU - u, IpU - u h ) + A h {u - u h , l!fu - u h ) 


<\\\IpU~ 


*III / p“-«II 


Therefore, 

and 


\IpU-u h \l < \\\I”u-u\\ 


II“ - “hill* ^ 111“ - IpU\\ 


Then, using Lemma 3.5, we get 

ML, III <" tmin (s—l.p) ||..|| 

111“ ~ “hill* ~ h ||“||H«(fi,ri)- 


(13) 


For the L 2 error estimate, we consider the following adjoint problem: find £ 
such that 


J — A£ = u — Uh, in J7, 

\ dnC = -«C + /3A r C + (“ - Uh), onTi, 

As u—Uh E L 2 (Ll, Ti), using Theorem [3j] yields an unique £ E H 2 (Q,T 1 ) 
satisfiying the following stability estimate 

IICII.ff 2 (Q,r 1 ) < 11“ — “h||L2(n, ri )- 


10 




(14) 


Using Lemma | .3.5 1 with p = 1, we get 

IIIC - /kill* ~ ^llell// 2 (n,ri) ;$ - «/i ||L 2 (o,ri)- 

Since Ah(-, •) defined in Q is symmetric, it is easy to see that it holds 
AhixX) = (u - «7i,x)L 2 (n) + {u ~ u h ,x)L2(Ti) Vxei / 2 (%). (15) 

Next, choosing x = u — in © and employing © together with ([9]) , 
we find 

ll 11 — u h\\L 2 (n,n) = Ah(u — Uh,C ) 

= Ah{u-u h , c-4c) 


The thesis follows using (13) and (14). 


□ 


4 The parabolic problem and its fully-discretization 


In this section we employ the results obtained in the previous section to 
present and analyze a DG space semi-discretization combined with an back¬ 
ward Euler time advancing scheme for solving the following parabolic prob¬ 
lem: 


dtu = A u + /, in 0, 0 <t <T, 

d n u = —au + /3Aru — Xdtu + g, on Ti, 0 < t < T, 

periodic boundary conditions, on T 2 , 0 < t < T, 

. u \ t= o = u 0 , in ft, 


(16) 


where T > 0, a, (3, A are positive constants and /, g, uq are (regular enough) 
given data. The weak formulation of (16) reads: for any t E (0, T\, find u 
such that: 


( {dtu, n) L 2 (n) + A (d t u, v) L 2 {Tl) + a(u , v) 
\«|t=o = ^ 0 ) 

for any v E i/ 1 (G,T i). 


Ui v )l 2 (Q.) + (9i v )l 2 ( ri)> 


(17) 


It is possible to prove the following result dealing with the existence and 
(higher) regularity of the weak solution of (16). 
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Theorem 4.1. If uq G H 2 {Pl, Ti), / G ^(0, T; L 2 (0)) and3 G ^(0, T; L 2 (ri)) 
and the following compatibility conditions holds 


1. m := A n 0 + /(0, •) G L 2 (Q), 

2. ni| ri := /3A r u 0 - d n u 0 - au 0 + g( 0, •) G T 2 (ri), 


then problem (16) admits a unique solution u with 

u G C([0, T]; H 2 (Q, Ti)) D ^([0, T]; L 2 (ft, r x )) n ^(0, T; H&Sl, Id)). 

Moreover, if uq G ^(fijld), g G iT^O, T; H 2m ~ 2k ~ 2 (n)) and J# G 
dd 1 (0, T; df 2m_2fe_2 (ri)), /or k = 0,..., m— 1 and the following higher order 
compatibility conditions hold for k = 1 ,,m 

3. u[ k) := Ai4*" 1} + |^/(0, •) G L 2 (Q) 

4 . := f3A r u^~^ - d„uf _1) - ■) G L 2 (Ti), 


where we set n^ 0) := u\ and = tt X | ri , tden it holds for k = 0 ,..., m — 1 

^ g c([o,T] ; i/ 2 m - 2 fe (o,ri))nc 1 (o,T;id 2 m - 2 fc - 2 (o,ri)) 

n ^(O.T;^"*- 2 *- 1 ^,?!)). ( 18 ) 

Proof. See Appendix A. □ 

Employing the DG notations introduced in Section |3.1[ the space semi¬ 


discretization of problem (16) becomes: find Uh G C°(0,T;V p (Th)) such 
that, for any t G (0, T], 


I (9 t u h ,v h ) L 2 (Q) + A {d t u h ,v h ) L 2 ( Tl) +A h {u h ,v h ) = (f,v h ) L 2 (n) + {g,v h ) L 2 {Vl) , 
= u h0, 

(19) 

for any Vh G V p {Th ), where Uho G V p (Th) is the L 2 -projection of no into 

V p (T h ). 

The following result shows the existence of a unique solution Uh of problem 

((191- 


Theorem 4.2. 

tion. 


The semi-discrete problem (19) admits a unique local solu- 
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Proof. As the proof is standard, we only sketch it. Let {(fj}j =1 be an orthog¬ 
onal basis of V p (Th )• The semi-discrete problem (19) is equivalent to solve, 
for any t E (0,T], the following system of ordinary differential equations 


j ( dt.Uh , 4>j)L 2 (vt) + dtUh , ^i)L 2 (ri) + -Ahiuh, <j)j) — + (<7> 0i)z, 2 (ri)> 

«/i|t=0 = u h0, 

„ ( 20 ) 


1 


for j = 1, Setting Uh = YliLi Q (£)<&> (20) can be equivalently written 

as 

( Mc(t) + Ac(t) = F (t), 

\c(0) 


= c°, 


( 21 ) 


where c(f) = (qO))^^, c° = (c?)^^ with and, for 

i,3 = 1, -jAr, 

Ajj = 4>j), Mij = (fa, 4>j)L 2 {n) + A(&, ^i)i/ 2 (ri)j 


T) — (/, 0i)l 2 (q) + ( 5 , <fo)z, 2 (ri)- 

Since the matrix M is positive definite and F(f) E L 2 (0,T;M Ar ) invoking 
the well known Picard-Lindelof theorem yields the existence and unique¬ 
ness of a local solution c E iL 1 (0, TV; M), i.e. Uh E H l (0, TV; V p (Th) C 
C'dO, TV]; V p (Th)) with T/v E (0,T], □ 


The next result shows the stability of the semi-discrete solution of (19). 


Lemma 4.3. Let Uh be the solution of (19). Then it holds 


KCOII 


£5(n,ro 


+ 


:dt< 


H' Ufe °llL 2 (n,ri) + / dl/lli 2 (n) + llslli^ri))^- ( 22 ) 


Proof. Choosing Vh = Uh in (19) and using (10) we get 


1 d 


2 dt " Uh " L l( n ’ r i) 


+ 


z < 
* 


(ll/llL 2 (n) + II^IIl 2 ^!)) \\ u h\\L 2 (n, re¬ 


using the arithmetic-geometric inequality and the Poincare-Friedrichs in¬ 
equality for functions in the broken Sobolev space iL 1 (7h), i.e., 


IklUw S(|o/.l^, (r , ) + IIMII 2 tl(£iiU£ r 1) ) 1/2 v h eH\%) 
IMU»(rO Stkl^^r^ + IIMII^r,) 1 ^ v h £H\T k ) 
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cf. 0 ], we obtain 

^IWli^ri) + IKIII* ~ ll/lll 2 (n) + ll^llz, 2 ( ri )- ( 23 ) 

The thesis follows integrating between 0 and T and noting that 

H Uw) IIl 2 (n,ri) ~ ll u °lli 2 (n,ri) because u ho is the L 2 -projection of u 0 into 

VP(T h ) X . ’ H 

Finally, we consider the fully discretization of problem ( |17[ ) by resorting 
to the Implicit Euler method with time-step At > 0. Let t k = kAt, 0 < 
k < K, with K = T/Af, and denote by u k ,k > 0,the approximation of 
u h (t k ). The fully-discrete problem reads as follows: given u° h = Uho , find 
u^ +1 G V p (Th), 0 < k < K — 1 , such that 

+ +A h (u k h + \v h ) (24) 

L 2 (fi) V / L 2 (ri) 

= (/(4+l), v h)L 2 (n) + (<?(4+l), ^)L 2 ( ri ) 

for all € V p (Th)- 



5 Stability and error estimates 


This section is devoted to show that the solution of problem (24) converges 
with optimal rate to the continuous solution of (16). We first prove the 
following stability result. 


Lemma 5.1. Let f k = f(t k ) and g k = g(t k ), k = 1, Then it holds 


k _ 


K 


II u h lll 2 (n,ri) + 


k=1 


K 


< 


k=1 


^oIIlUqA) + ^2 ( H^ / ' lli 2 (n) + lls^Hi^ro) • (25) 


Proof. We choose Vh = rt^ +1 in (24). Using (10), the identity 
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and the Cauchy-Schwarz inequality, we obtain 


U,k+1\\2 

I U h ll : L 2 (n, ri ) 


< A t 


fe|l L^(n,r 1 ) + II u h +1 II (n, ri ) + ^lll'“h +1 |ll* 

k+1 \\L 2 (n)\\ u h +1 \\L 2 (n) + ll5' /l+1 ||L2( ri )||n^ +1 || L 2( ri ) 


Employing Young’s inequality, Poincare-Friedrichs’ inequality and summing 
over k we get the thesis. □ 

We next state the main result of this section. 

Theorem 5.2. Let u G C([0, T]; H s (Ll, ri))n^ 1 (0, T; L 2 (fl, Pi)), s > 2, be 
the solution of (11) and let Uh be the solution of (24 )■ Ifdtu 6 L 2 (0, T ; H s (Ll, Ti)), 
6 i 2 (0,T; P(0,ri)) and satisfies 


_, 0 || <" i,min(s,p+l) I., II 

\ u 0 ~ ^11^2(0, ri ) ~ h IPoll H»(T h )’ 


(26) 


then 


L , K _ ,,X ||2 

\ u u h IlL^n.ri) 


<L 2 min(s,p+l) ( 11X112 1 II,. 112 

Up llx*(Q,ri) + ll n o||//s(n,ri) 


+ / \\dtu(t) ||//s(n,ri) ^ 


+ At 2 


lUt M (^)llL|(n,ri) 


and 


X 


At E 


11 111 2 

\\ u ~ u h\L 


k =1 


X 


,fc 112 


iHj[(n,ri) 


< tfrmn(s-l,p) I At ^ [ 

' fc=l 

+ ^ 2 |l u o|lHj(n,ri) + ll^' u ( i )llx*(n,ri) ^ 
lUt IIL|(f 2 ,ri) 


T 

2 / ||o2„ 


+ At 

where u k = u{tjf), k = 1,..., A". 

Proof. We first define the elliptic projection P : H 2 (Ll,Ti) —> V p (Th ) as 

A h (Pw - w, v h ) = 0 \/v h e V p (Th), (27) 
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where Ah{-, ■) is defined as in ([d]). We note (see Theorem |3.6[ ) that P satisfies 
the bound 


II Pw - w|| £ a (n)ri) + h\\\Pw - will,, < /i min(s ’ p+1 )||w||^ (n)ri) , 


(28) 


for all w E H S (Q, Ti), s > 2. We next write u k —u k = (u k —Pu k )+(Pu k —u k ) 
and start to focus on the second term. Considering problem (19) at time 
t k+ 1 , we easily get 


fPu k+l _ Pu k 
{ - At - 


) + A 

/ l 2 ( n) 


(Pu k+ 1 - Pu k 




At 


,v h | +A h (Pu 
l 2 (E) 


fc+i 


,Vh) 


— {f(tk),Vh)L 2 (Cl) + (g{tk), v h)L 2 (r!) — (E k+1 , Vh)L 2 (Q) ~ \{E k+1 , V^l 2 ^) > 

(29) 


for all Vh E V p (Th), where 


E k+1 = d t u(t k+1 ) - ^ (Pu k+1 - Pu k ). 


Subtracting (24) from (29), we get that e k = Pu k — u k satisfies 
: i Vh ] + A 


k+1 


7> " ' u “ft 
„k+l a k 


At 




L 2 (fl) 


At 

pk+l 


,Vh 


+ Ah{efr +1 , Vh) 


P(ri) 


= -(E k+ \v h ) L 2 {Q) - A (E k+1 ,v h ) L2{ri) , 


for all Vh E V p (Th )• Then, reasoning as in the proof of Lemma |5.1| , we 
obtain 


K 


II e h llz,2(o, ri ) + ^ ^2 lll e ftlL ~ H e ftllL2(r2,ri) + ^y^ll^llr^rn.riV ( 30 ) 


K 


k =1 


ft=l 


We bound the first term on the right-hand side of (30) using (26) and (28): 


IkftlL^n) — \\P u o ~ ^ftolL^n^i) 

< \\Puo - «o||L 2 (fi, ri ) + ll^o - «fto||L 2 (n, ri ) 

<*“»(«, P+I)||« 0 || w . (7i) , (31) 
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In order to bound the second term on the right-hand side of (30) we observe 
that it holds: 


E k+1 = 


Q ^ _ U k+1 + ( U k+1 - Pu k+1 ) ~ ( U k ~ Pu k ) 


At 


At 


‘ (t — tk) d 2 u(t) dt + — / dt(u(t) — Pu(t)) dt 

Jtk 


At 


tk 


where we employed Taylor’s formula. Therefore, employing the commuta¬ 
tion of the operators P and dt, we have 


| pfc+l112 ^ 

1 ~ At 


1 

+ A t 


rtk +1 

Ith 

rtk+i 


(t — tk) d 2 u(t) dt 


LPn,ri) 




( dtu(t ) — Pdtu(t))dt 


l\(P,v i) 


Using the Cauchy-Schwarz inequality we get 

C t k +1 


tk 


(t — tk)dfu(t) dt 


k+1 


< f / (t - t k ) dt 

>t k 

fi-k + l 


Ll(n,n) 

/ ftk+1 


r l k +1 \ 

J ll^t “Wllij^ri) 


1/2 


< 

r^i 


A * 3/2 (£‘ + ‘||a 2 »(t)|| 


L 2 (n ,ri) 


\ 1/2 

dt) 


Hence, 

1 

At 


Pk+1 2 

(t — tk)d^u{t) dt 


L l(n, ro 


Pk+1 

<At 2 WdMVWll^) dt - 

•ttk 
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Employing dtu € L 2 (0,T; H s {Th )), s > 2, and (28), we obtain 

rtk+i 


rtk+i 

/ ( dtu(t ) — Pdtu(t))dt 

Jt k 


Ll(n,r i) 

/ /-tfc+i \V 2 / /•4+i \ V 2 

“ ^ d *J (j t W d t u (t) ~ Pd Mt)\\ 2 Li { n, Tl ) dt 

, f /‘ ife + 1 \ 1//2 

< At 1/2 II d t u(t) - Pd t u(t) 11^2.0 r, ■> dt 

tfc+i \ H 2 


< At 1 / 2 /i min ( s,p+1 ) 




Il^«(t)||fl-s(f2 ) n) ^ 


Hence, 


1 

At 


ftk+l 


' tk 


( dtu(t ) — Pdtu(t))dt 


l\{ n,ri) 


f*fc+i 


< fc 2P+1 \\g,u(t)f H . (a Tl) dt. (32) 


Finally, summing over k we get 
K 


^ En®' 


fc=X 




(33) 


< At 2 J* l|8?..(t)lli 5(n , ri) dt + /! 2 “"<*-p+ ‘) J f |lft»(i)|l».( I! ,r 1 ) dt, 

a k 


which concludes the bound for e£. Finally, the thesis follow employing the 
triangle inequality and the bounds (30)-(31) together with (28)-(33). □ 


6 Numerical experiments 

In this section we present some numerical results to validate our theoretical 


estimates. In the first two examples (cf Sections 6.1 and 6.2) we consider 


a test case with periodic boundary conditions and validate our theoretical 


error estimates. In the last example (cf Section 6.3) we show that our theo¬ 
retical results seem to hold in the case of more general boundary conditions, 


provided the exact solution of problem (16) is smooth enough. 
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6.1 Example 1 


We consider problem (16) on Q = (0, l) 2 and choose / and g such that 
u = e" 10t (l — cos(2-7tx)) cos(47ry) is the exact solution. 

We have tested our scheme on a sequence of uniformly refined structured 
triangular grids with meshsize h = \/2/2 f , t = 2,...,7. In those sets of 
numerical experiments we have measured the error e(T) = u(T) — Uh.(T ) at 
the final observation time T = 0.001 in the ||- ||z, 2 (r2) an d IHIz, 2 (ri) norms. We 
have also measured the quantity Me^W*) 1 / 2 , being e k = u k — u\ . 

In the first set of experiments we used piecewise linear elements (p = 1) 
and the following parameters: a = 10, At = 10~ 5 , A = 10, f3 = 5 a = 2. 
The computed errors and the corresponding computed convergence rates 
are reported in Table [l] We have repeated the same set of experiments 
employing piecewise quadratic elements (p = 2); the results are reported in 
Table [2j From the results shown in Table [l] and Table [2j it is clear that the 
expected convergence rates are obtained. 


h 

\\ e ( T ) U 2 (f 2 ) 

rate 

l|e(T)|| L 2 (ri) 

rate 

(AtEAjrt.T 2 

rate 

y/2/2 2 

1.836048e-01 

- 

1.908256e-01 

- 

2.281359e-01 


V2/2 3 

5.455936e-02 

1.75 

5.035380e-02 

1.92 

1.186343e-01 

0.94 

v/ 2 / 2 4 

1.451833e-02 

1.91 

1.278655e-02 

1.98 

5.939199e-02 

1.00 

\/ 2 / 2 5 

3.688202e-03 

1.98 

3.208881e-03 

1.99 

2.962468e-02 

1.00 

V2/2 6 

9.258142e-04 

1.99 

8.028862e-04 

2.00 

1.480150e-02 

1.00 

V2/2 7 

2.316573e-04 

2.00 

2.006754e-04 

2.00 

7.399580e-03 

1.00 


Table 1: Example 1. Computed errors, p = 1, cr = 10, At = 10 5 , T = 0.001, 
A = 10, P = 5, a = 2. 


6.2 Example 2 

In the second example, we explore the dependencies of the error on the time- 


step At. To this aim, we set / and g as in Section 6.1 In Table [3] we report 


the computed errors and convergence rates obtained with piecewise linear 
elements (p = 1) and the following parameters: k = 7, cr = 10, T = 0.1, 
A = 10, /3 = 5, a = 2, h = V2/2 1 and vary the time integration step At. 
The numerical results are in agreement with the theoretical estimate. 
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h 

\\ e (T) L 2 (Q) 

rate 

\\ e (T) L 2 (ri) 

rate 

(AtEfrtlll^llll ) 172 

rate 

y/2/2 2 

2.470397e-02 

- 

1.751588e-02 

- 

5.281897e-02 

- 

\/ 2 / 2 3 

3.027272e-03 

3.03 

2.232268e-03 

2.97 

1.405198e-02 

1.91 

v/ 2 / 2 4 

3.827204e-04 

2.98 

2.822643e-04 

2.98 

3.602372e-03 

1.96 

V2/2 5 

4.797615e-05 

3.00 

3.539247e-05 

3.00 

9.081101e-04 

1.99 

V2/2 6 

5.992844e-06 

3.00 

4.421683e-06 

3.00 

2.276766e-04 

2.00 

V2/2 7 

7.507474e-07 

3.00 

5.632338e-07 

2.97 

5.593742e-05 

2.02 


Table 2: Example 1. Computed errors, p = 2, a = 10, At = 10 5 , T = 0.001, 
A = 10, $ = 5 a = 2. 


At 

e ( r ) L 2 (H) 

rate 

l|e(T)|| i 2 (ri) 

rate 

0.1 x 2 U 

2.682138e-02 

- 

8.678953e-02 

- 

0.1 x 2 _1 

1.487984e-02 

0.85 

4.905898e-02 

0.82 

0.1 x 2 -2 

7.889826e-03 

0.92 

2.630006e-02 

0.90 

0.1 x 2 -3 

4.050365e-03 

0.96 

1.360794e-02 

0.95 

0.1 x 2 -4 

2.028095e-03 

1.00 

6.881036e-03 

0.98 

0.1 x 2 -5 

9.897726e-04 

1.03 

3.415646e-03 

1.01 

0.1 x 2 -6 

4.664660e-04 

1.08 

1.656678e-03 

1.04 


Table 3: Example 2. Computed errors, k = 7, p = 1, a = 10, T = 0.1, 
A = 10, /3 = 5 a = 2. 


6.3 Example 3 


Finally, we consider problem (16) on A! = (0, l ) 2 with homogeneous Dirichlet 
boundary conditions applied T 2 and on Tj. In this case we choose / and 
g such that u = t(l — cos(27rx)) cos( 7 n/) is the exact solution. In Table 
[4] we report the computed errors and computed convergence rates at the 
final time T = 0.1. Those results have been obtained with piecewise linear 
elements (p = 1 ) and with the following choice of parameters: a = 10, 
At = 0.001, A = 10, /3 = 5 a = 2. We have ran the same set of experiments 
employing piecewise quadratic elements (p = 2 ); the computed results are 
shown in Table [5j The results reported in Tableland Table [5] clearly confirm 
the theoretical rates of convergence even in the cases of Dirichlet boundary 
conditions instead of periodic ones, at least whenever the exact solution is 


sufficiently smooth (see Remark 3.2). 
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h 

\\ e (T) L 2 (Q) 

rate 

\\ e (T) L 2 (ri) 

rate 

(At£f=, lll^lll.t 1 ' 2 

rate 

a/2/2 2 

9.185918e-03 

- 

1.111234e-02 

- 

1.347859e-01 

- 

/2/2 3 

2.704819e-03 

1.76 

2.849404e-03 

1.96 

6.413467e-02 

1.07 

/2/2 4 

7.279868e-04 

1.89 

7.169369e-04 

1.99 

3.155837e-02 

1.02 

/2/2 5 

1.875124e-04 

1.96 

1.797070e-04 

2.00 

1.571196e-02 

1.01 

V2/2 6 

4.745622e-05 

1.98 

4.501545e-05 

2.00 

7.847606e-03 

1.00 

V2/2 7 

1.192746e-05 

1.99 

1.127502e-05 

2.00 

3.922783e-03 

1.00 


Table 4: Example 3. Computed errors, p = 1, a = 10, At = 0.001, T = 0.1, 
A = 10, p = 5 a = 2. 


h 

e (^) L 2 (Q) 

rate 

\\ e (T) L 2 (ri) 

rate 

(At£t. 1 ll| e i lll. z ) 1/2 

rate 

a/2/2 2 

1.239177e-03 

- 

1.607590e-03 

- 

2.589798e-02 

- 

a/2/2 3 

1.543449e-04 

3.01 

2.189412e-04 

2.88 

6.771702e-03 

1.93 

/2/2 4 

1.911957e-05 

3.01 

2.788057e-05 

2.97 

1.715537e-03 

1.98 

a/2/2 5 

2.386211e-06 

3.00 

3.496808e-06 

3.00 

4.307186e-04 

1.99 

/2/2 6 

2.990873e-07 

3.00 

4.364171e-07 

3.00 

1.079691e-04 

2.00 

/2/2 7 

3.777961e-08 

2.98 

5.420558e-08 

3.01 

2.621607e-05 

2.04 


Table 5: Example 3. Computed errors, p = 2, cr = 10, At = 0.001, T = 0.1, 
A = 10, p = 5 a = 2. 


A 


Proof of Theorem 


4.1 


Proof of Theorem \f.l\ As the proof follows is based on standard arguments 
(see, e.g., m Chapter 7.1]), we only sketch the main steps. 


1. Construction of the discrete space. Let {e,}j>i be an orthonormal basis 
of L 2 (Q) such that 


Ve* • Vz = A i eiZ Vz e i > 1, 


i.e., A i and e* are respectively the eigenvalues and eigenfunctions of the 
weak form of eigenvalue problem —Ae = Ae with homogeneous Neumann 
and periodic boundary conditions on Ti and T 2 , respectively. Reordering 
{e«}*>i such that Ai = 0, it is easy to see that there holds 


Vej • Ve^ = 0, for i / j 


and 


/ |Vej| 2 = Aj > 0, for i > 1. 

Jn 
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Let V n = span{ej : i = 1, n > 1, and let tig be the L 2 (fl)- projection 

of uq on V n . Since the domain is regular, the eigenfunctions e* belong to 

H 2 (Q). 


2. Finite-dimensional approximation of We introduce the following 

finite dimensional problem: find u n G H 1 ^ 0, T ; V n ) such that, for t G (0, T), 


(dtu n , + X(dtu n , z ) L 2 ( ri ) + a(u n , 2 ) — (/, z) L 2^ + ( 5 , 

_ „. 7 i 

u |t=o — u 0 ) 


(34) 

for all ^ G F n , In the sequel we prove that problem (34) admits a unique 
solution in iL 1 (0,T; V n ). We write 


u n {t) = Uj(t)ej. 
3 =1 


The problem (34) is equivalent to find u (t) 
such that, for each t G (0,T), 


(«i {t),.:,u n {t)) T 


I Mu(t) + Au(f) = F(t), 

\u(0) = (lt 0) l, ...,U 0 ,n) T , 

where, for i, j = 1, ,.,n, 

Mij = + \M Vl := Sij + A(ej, 

Aij = a(ei,ej), Fi = (f,ei) L 2^ + (g,ei) L 2( Fl ), tto,i = (tto, ej)L 2 (Q)- 

Since the matrix M ri is semi-positive dehnite, we see that M is positive 
definite. In addition, F(f) G L 2 ( 0,T;M n ) and A : M n — >• M n is Lipschitz 
continuous. Therefore, by standard existence theory of ordinary differential 
equations, there exists a unique solution u(f) for a.e. 0 < t < T. 


3. Energy estimates. Taking 2 
inequality, we obtain 


u n in (34) and using the Cauchy-Schwarz 


d_ 

dt 


u 


1)) 


+ liv^f i2(n) 

< II 

11 


+ a\\u" ||| 2 ( ri ) + /3\\ Vrw' 1 j|^ / 2 ( ri ) 

u " Wiling) + ll/lli 2 (o) + llfflli 2 ^) 


(35) 
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for a.e. t £ [0,T]. Using the differential form of the Gronwall’s inequality, 
data regularity and Lemma A.l we obtain 

< 11 «0 11 L2 (n, ri ) +11 /11 L 2 (0,T;Z,2 (fi)) +11 5 11 i 2 ( 0 ,T;L 2 (Fi)) — ^' 


Integrating (35) in [0 ,T] and employing the above inequality together with 
data regularity and Lemma [A .1 1 we get 


l un |lL 2 (0,T;//i(a,ri)) ~ ll“o|| L 2( n)ri ) + ||/|| L 2( 0)T;i 2(n)) + ||5llL 2 (0,T ; L 2 ( ri )) - 


< c. 


On the other hand, taking 2 = dtu n in (34), integrating in t and using the 
Cauchy-Schwarz inequality, we obtain, for every r £ (0, T], 

\ r IMlWro + I||V«“(r)||| 1(n) + % l|u"(r)|| 2 i3(ri) + §||Vr»"(r)& (ri 


1 ,, — „ ,, o Oi ,, n 11 2 , P ] 


2 
77-11 2 


< - 2 \\V<\\Un) + ^IKIIl 2 ^) + 5l|Vr<|ll 2(n) 

2 J 0 l|/|l * 2 («) + 2A J 0 ll5ll ^ 2 (ri)’ 


2 

1 

+ 77 


where the right-hand side of the above inequality can be bounded using 
Lemma | A. 1| and data regularity. 

Moreover, differentiating (34) with respect to t and setting u n := dtu n 
we get for any t £ [0, T] 

(dtu n , z)l 2 ( ci) + Hdtu n , z)l 2 (ri) + a{u n , z) = (dtf, z)l 2 (q.) + {dt9, z) L 2 ( ri ), 

(36) 

for all 2 £ V ". Testing (36) with z = u n , it is easy to show that it holds 

pt 


rt 


\\®t u IIz-2(n,ri) + / II ®t u ( s )ll_H'Ho,ri) ~ / ll^/( s )llL 2 (n) ds 


0 


+ / IIStsMIl L 2 (ri) ds + ||^(o)||i i(nri 


)• 


(37) 


Taking t = 0 in (34), testing with 2 = dtu n { 0), integrating by parts and 
employing the Cauchy-Schwarz inequality once more, we obtain 

ll^ n (0)lli 2 ( n, ri) IKWII^iq) + 11/(0, 01112(0) + Il9(0,-)lli 2 (r l)5 

whose right-hand side can be bounded by resorting to compatibility condi¬ 
tions, Lemma A.l and data regularity assumptions. 
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Hence, collecting all the above results, we get 

u n € c([ o, T]-Hl(n, r,)) n c\o, t- l\{^ ro) n ^(o, t- h\$i, id)). 

4. Existence of the solution u. Resorting to subsequences {u mi }ff zl of 
{ u m}m=i , passing to the limit for m —> oo and using standard arguments it 
is possible to prove that there exists a solution u to problem © with 

u £ C([0, Id)) n C' 1 (0, T; n)) D H\ 0, T; Pi)). 


5. Uniqueness of the weak solution. Let u\ and U 2 be two solutions of weak 
problem © and set w = u\ — U 2 - By definition, taking z = w, we get from 

© 

fff (ll' u/ llz,2(f2 ) n)) + ll Vu ’lli 2 (n) + a IMIi 2 (ri) + / 3 ||VrHll 2 (ri) = °> 

that implies w = 0, or u\ = U 2 for a.e. 0 <t<T. 

6. Improved regularity. Rewriting © as 

a(u,v) = (f,v) L 2 (n) + (g,v) L 2 { n), 


where f = f — dtu £ L 2 (0, T, L 2 (fi)) and g = g — dtu £ L 2 (0, T, L 2 (Id)). 
Employing Theorem 3.1 we get u(t) £ id 2 (f2, Ti) for a.e. 0 < i < T. 


6. Higher regularity. We prove (18) by induction. From the above discus¬ 
sion the result holds true for m = 1. Assume now the validity of (18) for 
some m > 1, together with the associated higher order compatibility and 
regularity conditions. Differentiating (16) with respect to t, it is immediate 
to verify that u = dtu verifies 


d t u = A u + f, 

d n u = —au + f3Aru — Xdtu + g, 
periodic boundary conditions, 

„ «|t=o = «o, 


in D, 0 < t < T, 
on R, 0 <t <T, 
on T 2 , 0 < t < T, 
in Q, 


(38) 


where / = d t f, g = d t g, u 0 = /(0, ■) + Au 0 in D and u 0 |r = /3A r u 0 - 
d n u o — crao + <?(0, •) on T. Since the pair (f,g) satisfies the higher order 
compatibility conditions for k = 1 ,m then the pair (/,<?) satisfies the 
same type of compatibility conditions for k = 1,..., m — 1. Hence, it follows 
for k = 0 ,..., m — 1 

^ £ C([0,T];lL 2m - 2fc (D,r 1 ))nC 1 ([0,T];F 2m - 2fc - 2 (L!,r 1 )) 

nff'fO ,T;iTj* m - 2fc “ 1 (fi,ri)) (39) 
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which immediately implies the validity of (18) for k = 0,... ,m. 


□ 


The following result has been proof in m Lemmas 4.4 and 4.5]. 

Lemma A.l. Let z G Z = {z G H 2 (Ll) | d n z = 0 on Ti}. If z n is the 
L 2 (Q) -projection of z on V n , then 

|| z n — -s||i^i(r 2 ;r 1 ) ~^ 0 when n — > oo. (40) 

Let Voo = \Jff =1 V n . Moreover, Z and are dense in i). 
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